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SUMMARY 


The experimental wave drags of bodies and wing-body combinations 
over a wide range of Mach numbers are compared with the computed drags 
utilizing a 24 -term Fourier series application of the supersonic area 
rule and with the results of equivalent -body tests. 

The results indicate that the equivalent -body technique provides a 
good method for predicting the wave drag of certain wing-body combina- 
tions at and below a Mach number of 1. At Mach numbers greater than 1, 
the equivalent -body wave drags can be misleading* The wave drags com- 
puted using the supersonic area rule are shown to be in best agreement 
with the experimental results for configurations employing the thinnest 
wings. The wave drags for the bodies of revolution presented in this 
report are predicted to a greater degree of accuracy by using the frontal 
projections of oblique areas then by using normal areas. A rapid method 
of computing wing area distributions and area-distribution slopes is given 
in an appendix. 


INTRODUCTION 


The area rule, first advanced by Whitcomb in reference 1, has con- 
siderably altered the methods for predicting wave drag of wing-body com- 
binations* Studies leading to the discovery of the area rule showed that 
interference drag between wing and body components could be very large. 
Therefore, estimation of drag by component buildup without somehow evalu- 
ating the interference drag could give misleading answers. However, in 
consequence of the transonic area rule, a valuable tool was ma.de avail- 
able to the designers in assessing the transonic drag. This was the 
equivalent-body concept, which states that at transonic speeds the pres- 
sure drag of the airplane is the same as that for a body of revolution 

^Supersedes declassified NACA RM L56D11 by Robert L. Nelson and 
Clement J. Welsh, 1957- 
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having the same longitudinal distribution of cross-sectional area. As a 
result, the drag of the configuration is obtained by either estimating 
or experimentally determining the equivalent -body drag. Experimental 
checks for airplane configurations presented in reference 2 generally 
support this concept in the transonic speed range. 

The supersonic area rule, given by Jones in reference 3, provided a 
powerful method for calculating the wave drag at supersonic speeds. In 
references k and 5 the mechanics of the drag calculations were discussed 
together with a number of comparisons of calculated and experimental 
drags generally at low supersonic speeds. Jones pointed out in refer- 
ence 3 that the method could be expected to give good results for thin 
wings mounted on vertically symmetrical bodies. Later, Lomax in refer- 
ence 6 gave the complete linearized theory expressions for the drag. 

The added terms in Lomax 1 s result represented the limitation pointed out 
by Jones. 

The purpose of the present paper is to provide a better feel for 
the range of applicability of both the transonic and supersonic area 
rules. For the transonic area rule, this is done by making additional 
comparisons between equivalent -body and wing-body experiments. For the 
supersonic area rule, comparisons are made of calculated and experimental 
results for both body and wing-body combinations over a wider range of Mach 
numbers than heretofore made. The supersonic-area-rule calculations were 
made by using a 2k -term Fourier series expression for the slope of the area 
distribution. 


SYMBOLS 


frontal projection of the area cut by a Mach plane or wing 
aspect ratio 


a n = — 

n Jt , 


C D 

C P 

c 

C 0 


r ( §a _ 

p f\ 

/ 0 

2 q] 


drag coefficient, 
pressure coefficient 
wing local chord 

root chord of particular pointed wing tip 


wing root chord 
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c t 

D 

d 

f 

f(v) 


wing tip chord 
drag 

maximum body diameter, 2r m 
resultant pressure force 

wing-thickness-distribution function, - - ■ 

Z max 


G(K,v 0 ) 

H(K,v 0 ) 

K 


l 

h 

ij d 
M 


wing-area-distribution function 
wing-area-distribution slope function 


mfl + — CQS - for the left-wing panel; m(l 
V tan A / \ 


tan A 

the right-wing panel 


|3 cos 0 
tan A 


length of configuration 

total length of area distribution 


body fineness ratio 
Mach number 


m = ^ ^ + tan A 
£ (1 - A) 

n integer 

q dynamic pressure 

r local body radius 

r m maximum body radius 


S 

*b 


reference area 
body frontal area 

wing exposed area 


^ for 


wing total plan-form area 



h 


s 

wing semispan 


s o 

semispan of particular pointed-tip wing 

t/c 

wing thickness ratio, 

2z max 

c 


Cartesian coordinates 


x o 

point of intersection 

of Mach plane with the x-axis 

X* 

x-coordinate measured 

from wing leading edge 

X 1 

dummy variable 



z local maximum wing ordinate 


p = /m 2 - 1 



0 = cos’^( 2 — - 1 

V l t 


A 


A 


wing leading-edge sweepback angle , except as indicated other- 
wise by subscripts c/4, c/2, and 3c/4 for sweep at other 

chord lines 

c+ 

wing taper ratio, — 

c r 

Mach angle, sin~^ i 


v 


c 


Vq value of v at root of particular pointed-tip wing 

v r value of v at root of actual wing 


1 

] 

C 

c 

c 


0 


angle between z-axis and line of intersection of Mach plane 
with the y,z plane 


REVIEW OF THE BASIC THEORY 


From reference 6, the equation for the wave drag of any system of 
bodies or wings and bodies can be written as 
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D 

q 



dX n 


d 2 A(x,9) 


dx 


_P_ df(x,9) 
2q dx 


d 2 A(x-| ,9) 



2q dxj 


log (x 



( 1 ) 


The equation is subject to the usual limitations of the linearized theory. 


Before discussing the terms in the drag equation, it is well to review 
the definition of Mach planes. The physical significance of equation (l) 
is understood if the configuration is cut by Mach planes. Mach planes are 
easily visualized by considering a Mach cone originating at a point on the 
x-axis which is alined with the remote relative wind. A Mach plane is 
simply a plane tangent to the Mach cone and at an angle of roll, 0 about 
the x-axis measured from the y-axis. By moving the vertex of the Mach 
cone along the x-axis, a series of parallel Mach planes will cut the con- 
figuration for a fixed roll angle 0. 


In the drag equation the term A(x,0) represents the frontal pro- 
jection of the oblique area cut by a particular Mach plane, whereas f(x,0) 
represents the net force normal to the stream direction on this section in 
the 0 direction. These relationships are illustrated in figure 1 for 
angle of roll 0 of the Mach plane of 0° and 90°. By neglecting the 

term ^ the equation reduces to the supersonic-area-rule formula 

given by Jones in reference J. Evaluation of f(x,0) requires the pres- 
sure distribution on the conf iguration which when integrated over the 
configuration gives the drag directly. As a result, large values of 

^ impose a limitation on the supersonic area rule, even within 

the framework of the linearized theory. 


It is not the purpose of the present paper to evaluate the drag of 
configurations including the effect of the pressure term but to evaluate 
the drag of configurations using the supersonic-area-rule formula of 
Jones. The influence of the pressure term was evaluated for one simple 
case. 


Equation (l) can be written in coefficient form as 



( 2 ) 
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where 


Cj)(Q ) — ~ 


nL rl 


t- l A 2 


2rtS 


0 u 0 


4 - -J- 4^- ] [ ^ A 9 - — — jlog (x - x^)dx dx^ (3) 

y dx 2 2q dxJl dXi 2 2q dx ± ) V ±; 1 


The quantity Cj)(0) is most readily determined by solving the inte- 
gral for C j~)( 0 ) through a Fourier sine series expression for ^ fol- 
lowing the method of reference 7 if 


91 = ““'H 2 H ■ x ) 

(# - 1 ^ sin ■* •* 


then C-p(0) can be written as 


c d( 9) = 17^2na n 


4s 

For the computations of this paper, only 24 terms were used in the 


Fourier sine series expression for 


dx 


Thus. 


n=24 


c d( 0 ) = j^T Y, *** 


n=l 


BODY DRAG RESULTS 


For bodies of revolution, the calculation of the drag is simplified 
to some extent because the area distributions are identical for all roll 
angles. However, except for high-fineness-ratio bodies, it is not pos- 
sible to assume that the frontal projection of the oblique area cut by 
the Mach plane is the same as the normal area. Figure 2 shows an example 
of this for two parabolic bodies of revolution having different fineness 
ratios and shapes. The area-distribution curve slopes were calc ula ted 
from the expression 

dA _ _2_ f x u (x - xp)dx 

dx - (x - xq ) 2 
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The derivation of this expression is given in appendix A. It has also 
been assumed for the calculations (and all succeeding body calculations) 
that a cylinder can be added at the base of the body without altering 
the drag. If this were not done, the solution would require the flow to 
fill the area behind the base which would exceed the limitations of the 
linearized theory. Figure 2 shows large changes in the peak slope over 
the afterbody of the fineness-ratio-6.04 configuration; these changes 
would lead to a significant drag variation with Mach number. 

The evaluation of the slope of the oblique area distributions is 
extremely difficult except for simple bodies. There naturally arises 
the question as to whether this is worth while if the pressure term is 
ignored. 

As derived in appendix B, the local force acting on the oblique 
area of a body of revolution is 




r dC^ 

The only assumption made in derivation of — is that — is 

1 dx 

constant over the oblique area. This is a reasonable assumption except 
for bodies having discontinuities, and high local slopes. Then, 


dA _ P f _ dA + ££p 

dx 2 q dx 2 dx 


Thus, the error in the drag introduced by ignoring the pressure term is 

dCp 

dependent on the pressure gradient — =-. 

dx 


It would be expected that the drag for a conical nose with an 
attached shock wave over which the pressure is constant at zero angle of 
attack would be least affected by the pressure term. (The pressure term 
takes on a value only near the juncture with the cylinder; however, the 
pressure term was not evaluated in this region.) Figure 3 presents a 
comparison of the drag of various cones calculated with the supersonic- 
area-rule formula with the exact theory drag of reference 8. The lowest 
Mach number of the comparison corresponds to the lowest Mach number for 
entirely supersonic flow on the cone as calculated with the exact theory. 
The highest Mach number of the comparison was arbitrarily taken as that 
at which the slope of the Mach line equaled one-half the slope of the 
nose. The agreement between the two theories is remarkable, within 
5 percent except for a few points. 



8 
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A better comparison may be made by plotting C D (z/r m ) against 
P r m/^-> q uant ity which defines the frontal projection of the oblique 

area distribution. This has been done in figure 4 to give drag in dimen- 
sionless or collapsed form. The drags from exact theory (5° half angle 
cone was chosen as representative), slender -body theory (ref. 9), and 
supersonic-area -rule theory are shown. The comparison shows the great 
improvement of the area-rule theory over the slender— body theory at 
values of |3r m /Z greater than 0.2, and the good agreement of the area- 

rule result with the exact theory to (3r m /Z of about 0.7. At higher 
values of $v m Jl the area-rule theory is in error, possibly first because 
the pressure term is neglected but finally, near |3r m /Z = 1, because the 
assumptions of the linearized theory are violated. At Pr m /Z = 1, the 

Mach line lies on the cone surface, which corresponds to the realm of 
hypersonic flows. (See ref. 10.) 

For a body with curvature, for example, a nose of parabolic profile, 
the pressure over the nose is variable, and the influence of the pressure 
term may be significant. Figure 5 presents the drag for noses of parabolic 
profile in collapsed form. Here the supersonic-area-rule theory is an 
improvement over slender-body theory but in only partial agreement with 
the more exact second-order theory of reference 11. Inclusion of the 
pressure term, evaluated by using second-order pressure distributions, 
however, does give agreement with some of the second— order— theory results. 
Since the second— order— theory drags do not collapse into one curve, agree- 
ment should be expected only with those points for which the pressure dis- 
tribution used in evaluating the pressure term apply. However, this was 
not the case. For example, the pressure distribution used for the pres- 
sure term calculation at Pr m /Z = 0.3 corresponds to the flagged symbol. 

For the parabolic noses, both the area-rule theory and the area-rule 
theory plus the pressure correction cannot be expected to apply near and 
above Pr m /Z =0.5, where the slope of the Mach line equals the slope of 
the nose tip. 

Figure 6 presents a comparison of the pressure drag from supersonic- 
area-rule theory with experiment and slender-body theory for a family of 
parabolic bodies of revolution. The experimental drags were taken from 
references 12 and 13» In determining the experimental pressure drags, 
the friction drag was assumed turbulent and evaluated by using the sub- 
sonic drag level and the results of reference 14 for the effects of Mach 
number and Reynolds number; the fin pressure drags were assumed identical 
and have been published in reference 15 ; and the base drags were small 
and were subtracted when available. The slender-body-theory drags were 
calculated using the curves of reference 9. 
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As would be expected, the comparisons show the increasing ability of 
both the area-rule theory and slender-body theory to predict the drag as 
the body fineness ratio is increased. In most cases, the area-rule theory 
offers a significant improvement over slender-body theory. The area-rule 
theory and slender-body theory are in agreement near M = 1, since at this 
Mach number the supersonic-area-rule theory reduces to slender-body theory. 

From these nose and complete-body comparisons that have been made, 
the following conclusion can be drawn. The area-rule drag of bodies can 
be predicted to a greater degree of accuracy by using the frontal projec- 
tion of oblique areas at a given Mach number than by using normal areas, 
if, at the Mach number under consideration, the limitations of the line- 
arized theory are not exceeded . This is illustrated by the comparison 
between the drag at a given Mach number and the drag near M = 1 especially 
for the low-fineness-ratio bodies. It is not to be inferred from the above 
statement that the supersonic -area-rule method is recommended for evalu- 
ating the drags of bodies of revolution. However, when the drags of wing- 
body combinations for which the body area distribution is needed are deter- 
mined, the oblique area distribution should be used if the body is of low 
fineness ratio or has low-fineness -ratio components. 


CALCULATION OF WING-BODY DRAG 


The difficulty in computing the wave drag of wing-body configurations 
can be considerably reduced if the configuration meets the following con- 
ditions: first, the body is of sufficiently high fineness ratio so that 

the change in body-area distribution with Mach number is small, and sec- 
ond, the wing is thin. These conditions imply also that the pressure 
term is negligible. Some feel for the body fineness ratios necessary for 
the above condition to be met can be obtained from the preceding section 
on bodies of revolution. The assumption of a thin wing allows the Mach 
plane intersecting the wing obliquely to be replaced by a plane perpen- 
dicular to the wing chord plane intersecting the wing plane along the 
same line as the Mach plane. Note that, at zero roll angle, the Mach 
plane is normal to the wing chord plane but is not normal to the wing 
chord plane at any other roll angle for a Mach number other than M = 1. 
Also the angle between the Mach plane and the normal to the wing chord 

plane is greatest and equal to tan -1 p at a roll angle of 90°. 

Appendix C presents a simple analytical method for evaluating wing 
area distributions and area-distribution-curve slopes. The curves neces- 
sary for evaluating these quantities (figs. l6 and 17) are applicable only 
to 65A series airfoils, but similar curves can be made up for other air- 
foil sections. 
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In order to get an idea of the applicability of the thin-wing assump- 
tion, a calculation has been made of the true area-distribut ion-curve - 
slope variation for 60^ delta wing having an NACA 65 AOO 6 airfoil section 
for a roll angle of 90° and a Mach number of 1.4l4. In order to simplify 
the calc ula tion, the wing was approximated by a sufficient number of 
ping ar— slope elements to define the airfoil section adequately* With 
this approx ima tion the Mach plane intersection with the wing surface was 
made up of straight lines. The expression for the frontal projection of 
the oblique area was then easily evaluated and differentiated to obtain 
the slope. The results of the calculation axe presented in figure 7* 
Although the slopes for the upper and lower half wings are significantly 
different, the total slope agrees almost exactly with the slope obtained 
by using the thin-wing assumption. On the basis of this result, it is 
felt that the thin-wing solution should be adequate for wings of present- 
day interest. 

For the wing-body combinations of this paper, an additional simpli- 
fication was allowed in the supersonic-area-rule wave-drag calculations. 
Since the tail fins mounted on the models were thin and relatively small 
(see ref. l 6 ), their drags were subtracted as tares. Then, since the 
bodies for all cases were of high fineness ratio (and identical), the 
body-area-distribution-curve slopes were considered independent of Mach 
number, and the changes in the area— distribution— slope curves with Mach 
number and roll angle were due entirely to the wings. As derived in 
appendix C, the area distribution for a given wing (m fixed) is depen- 
dent only on the value of P cos — . Thus, the area-distribution-slope 

tan A 

curves for the wing-body configuration are dependent only on the value 

of £ cos 9 . Then, from equation (2) and because of the symmetry of the 
tan A 

configuration. 


on/ 2 

C D (M) = § / C D ( 0 )d 0 

n J 0 

In order to obtain the wave drag of the configuration, a plot of 

against 0 is required. This can be computed if a plot of C D 

P cos - 9 - is given, since the angle 0 is known for fixed values 
tan A 

and P cos The configuration drag is simply the average drag 
tan A 

0=0 and i. 

For the wing— body calculations of this paper the bodies were identi- 
cal. The body-area-distribution-slope curves are shown in figure 2(b). 

The curve for M = 1.4l4 was chosen as representative for the Mach number 


against 

of —£-r 
tan A 

between 
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range of interest. The wing area distributions and area-distribution- 
curve slopes were obtained by using methods similar to that given in 

appendix C. In addition, a limiting value of £ -- ° - s 9 = 0.8 was set for 

tan A 

configurations having blunt leading-edge airfoils. (Above ft c . os 9 = 1 . 

\ tan A 7 

the Mach line lies behind the wing leading edge, and the linear theory is 
no longer valid for blunt airfoils 

/ 

An example of the wave drag calculation for the most extreme configu- 
ration investigated (60° delta wing, NACA 65 AOO 6 airfoil) is presented in 

figures 8 , 9> and 10. Figure 8 shows nondimensional plots of — against 

dx 0 

0 for various values of £ C0S Figure 9 shows the effect on Cr> of 

tan A u 

the number of terms in the series solution. Except at L- 92 s — ft = 0 and 0.8, 

tan A 

convergence was apparently obtained within 2k terms. Figure 10 shows the 

variation of the area distribution drag with P cos 8 -the variation of 

tan A 

area distribution drag with roll angle, and the variation of the config- 
uration drag with . ^ . 

tan A 


WING-BODY DRAG COMPARISONS 


Figure 11 presents some wave-drag comparisons for wing-body combina- 
tions. The experimental wing -body results were taken from references 16 
to 19* The wing -body wave drags were obtained in the following manner. 

The friction drags were assumed to be turbulent and were estimated by 
using the results of reference l4. Base drags and fin pressure drags were 
subtracted using the results of reference l 6 . The equivalent-body drags 
for a Mach number of 1 were obtained experimentally by using the helium- 
gun technique described in reference 2. These models had four scaled 
tail fins. The friction drag was assumed to be the subsonic drag level 
corrected at higher speeds for Reynolds number and Mach number by using 
the results of reference 14. Base-drag rise and fin-drag rise were not 
evaluated for the equivalent -body models. These quantities, however, 
should be small in the Mach number range where comparison is valid. The 
supersonic-area-rule-theory drags were evaluated by using the method of 
the preceding section. No attempt was made to evaluate the drag with 
the pressure term included. The drag coefficients presented in figure 11 
are based on total wing area. 

The inability of the supersonic-area-rule theory to predict the drag 
near M = 1 is evident for nearly all cases. However, the agreement at 
the higher Mach numbers between the theoretical drags and the experimental 
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wing-body drags is excellent and within the accuracy of evaluating the 
experimental wave drag, except for three conf igurations . Two of these 
configurations (figs. ll(c) and 11(g)) had 6-percent -thick wings which 
were the thickest wings investigated. The third configuration (fig. 11(h)) 

had a percent -thick airfoil but with fairly steep wedge components. 

For these configurations, a significant effect of the neglected pressure 
term may be possible. As a result, the drags calculated for configura- 
tions having wings of these thicknesses and sections should be viewed 
with caution. 


The comparisons in figure 11 show that the equivalent -body drags 
give a good approximation to the experimental wing-body drags up to a 
Mach number of 1, except for the two configurations having 6-percent- 
thick wings (figs. 11(c) and 11(g)). This result is in agreement with 
reference 2 which shows the validity of the transonic-area rule decreases 
with increasing wing-thickness ratio. At Mach numbers above 1, the agree- 
ment is variable but tends to be consistent with the flatness of the cor- 
responding theoretical curve. That is, as the theoretical drag variation 
with Mach number becomes smaller, the equivalent body gives a better 
approximation of the supersonic-drag level. This would be expected, since 
a flat theoretical curve indicates that the variation in area-distribution 
drag with roll angle or Mach number is small. Then the drag for the Mach 
number 1 or roll angle 90 ° area distribution (corresponding to the equiv- 
alent body) is representative of the configuration drag. 

Figure 12 shows the comparison between experimental-configuration 
drag and equivalent— body drag for an airplane configuration. The com- 
parison shows an extreme example, compared with the relatively good results 
of reference 2, of the inability of the equivalent body to predict the 
supersonic-drag level. The equivalent -body drag is approximately 40 per- 
cent low in spite of the low aspect ratio of the configuration. Apparently , 
the configuration tail surfaces cause the area distribution to change mark- 
edly at low supersonic speeds. Below M = 1, the equivalent body gives a 
fair representation of the configuration drag. The drag of the configu- 
ration minus the tail surfaces could probably be calculated to the degree 
of accuracy shown in figure 11. The influence of the tail surfaces, how- 
ever, may be difficult to evaluate. If the horizontal tail supports a 
load when the conf iguration is at zero lift, the influence of the pres- 
sure term may be significant. Although no supersonic-area-rule drag cal- 
culations were made for this airplane, reference 20 indicates that gen- 
erally good predictions of complete-airplane drag can be made.'*" 


^"Subsequent to the preparation of this paper, NACA RM A 56 IO 7 has 
been prepared at the Ames Laboratory and presents supersonic -area-rule 
calculations for a configuration similar to the one shown in figure 12 
but with small differences in area distribution in addition to the 
absence of a canopy . 
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EVALUATION OF COMPONENT AND INTEKFERENCE DRAGS 


As was shown in the preceding section, the supersonic area rule can 
be a useful tool in evaluating the supersonic drag of a wing-body con- 
figuration. In order to assess the efficiency of the combination as a 
whole, however, the effects of the combination on the component drags 
and the interference drag between the components must be known. The 
supersonic area rule provides a valuable method for evaluating these 
effects . 

The supersonic-area-rule equation can, of course, be used in evalu- 
ating the drags of individual wing and body components. This was done 
for a number of bodies in a previous section of this paper. The same can 
be done for isolated wings. An example of this is shown in figure 13 
where the drag of delta wings having 63 A series sections is plotted in 
collapsed form. The area-rule result is compared with a result obtained 
by the method of Beane (ref. 21). The two methods are just two forms of 
the same linearized wing theory. The agreement between the two methods 
is good. 

An example of the effect of the wing-body configuration on wing drag 
is presented in figure 14. The calculation is for the configuration having 
the closest agreement between the theoretical and experimental drags 
(fig. 11(b)). In this figure, the drag of the exposed-wing panels based 
on total and exposed wing areas is compared with the isolated wing drag. 
Separation of the wing panels gives approximately a 10-percent reduction 
In wing drag coefficient at Mach numbers above 1.3. As the Mach number 
approaches 1, this favorable effect disappears. This would be expected, 
for at M = 1 the area distributions of the exposed wing panels and the 
Isolated wing would be identical if the body were cylindrical. 

Figure 15 shows an evaluation of the interference drag for the same 
configuration. The sum of the calculated body and wing drags is compared 
with the calculated configuration drag. The curves show a favorable inter- 
ference effect at Mach numbers below M = 1.3- At Mach numbers above 1.3* 
interference drag is, for all practical purposes, zero. Thus for this 
configuration, at Mach numbers greater than 1.3* the only beneficial effect 
of combining the wing with the body comes from the separation of the wing 
panels . 


CONCLUSIONS 


An investigation has been made of abilities of the equivalent -body 
technique and a 24-term Fourier series application of the supersonic- 
area-rule method to predict wave drag at transonic and supersonic speeds. 
From the theoretical and experimental comparisons made, the following 
conclusions can be drawn: 
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1. The area-rule drag of the bodies of revolution presented in this 
report are predicted to a greater degree of accuracy by using the frontal 
projection of oblique areas at a given Mach number than by using normal 
areas . 

2 . The supersonic wave drag of slender-wing — body configurations can 
be predicted with the supersonic-area-rule formula. For the wing -body 
configurations investigated, the best agreement was obtained for the con- 
figurations employing the thinnest wings. 

3. The equivalent body technique provides a good method for predicting 
the wave drag of certain wing-body combinations at and below Mach number 1. 
At Mach numbers above 1, the equivalent body wave drags can be misleading. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 6, 1956. 


APPENDIX A 


AREA DISTRIBUTION SLOPE FOR BODIES OF REVOLUTION 
CUT BY OBLIQUE MACH PLANES 


The area distributions are identical for all roll angles. For 
simplicity a roll angle of 90 ° vill be used in the derivation. 



The frontal projection of the oblique area cut by the Mach plane (see 
sketch (a) ) is given by the equation: 


A 



y dz 


From the equation for the Mach plane, z is related to x by 


z = X ■" J S? 
P 

and 

dz = i dx 

P 


The equation relating y and x is given by 



/r 2 (x) - 


(x - xo)‘ 

1 ~ 


= T - (x - xq)‘ 



1 6 


Then, 


A = 2 


' u /p 2 R 2 (x) - (x - x 0 ) 2 dx 


the expression for A gives 


(A2) 




(* ~ xp) 

v/3* R 2 (x) - (x - Xq) 2 


dx 


(A3) 


where x z and x u are the roots of x = x Q - pR(x) and x = x Q + pR( x ), 

respectively. In the text the subscripts on x are omitted when it is 
not important to stress the functional nature of the symbol. 
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APPENDIX B 

THE NET FORCE ACTING ON THE OBLIQUE AREA OF A BODY OF 
REVOLUTION AT ZERO ANGLE OF ATTACK 


For a body of revolution at zero angle of attack, the net force is 
independent of roll angle. The derivation will be made for a roll angle 



The net force in the 9 direction (the z-direction for a roll angle of 
90°) (see sketch (b)) can be written as 


f 

q 



sin a dS 


Since dS = /dy 2 + dz^ and sin a = 




^dy^ + dz^ 


I ~/ c Cp a y 

The pressure coefficient at zero lift is a function of x only. 
The equation for y in terms of x is given by equation (l) of appen- 
dix A as follows : 


y = - (x - xq) 2 
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and 


dy 


P 2 R(x)^ - (x - xq) 

P \/3 2 R 2 (x) - (x - xq) 2 


dx 


Then, the net force can be written as 





lower surface 



C p[ p2R(x) H ~ (x ~ Xq) ] dx + 2 [* o c p[ p2R(x) ti ~ (x ~ *°j] 

^P 2 R 2 (x) - (x - xq ) 2 \/p 2 R 2 (x) - (x - Xq ) 2 

f = 2 fxu Cp[p 2 R(x);£ - (x - xoj] 
q P * l y/p 2 R 2 (x) - (x - xq ) 2 


Integrating the expression by parts gives 


£ 

<1 




yjffrPix) - (X - Xq ) 2 


dx 


jf £ j_s essentially constant between the limits of integrations, 

dx 

r u y s ¥(x) . (x . *0)" ax 

4 f> dx / » P 

Then, from equation for the frontal projection of the oblique area 
in appendix A (eq. (A2)) 


.p dC^ 

£ = - PA — £ 
q dx 
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APPENDIX C 

METHOD FOR DETERMINING WING-AREA DISTRIBUTION AND 
AREA-DISTRIBUTION-CURVE SLOPE 


This method assumes that the wing is thin and that the oblique Mach 
plane can be replaced by a plane perpendicular to the wing chord plane . 
The method also assumes the wing has straight leading and trailing edges 
and constant thickness ratio. 

The method is developed first for pointed-tip wings. Then, correc- 
tions are made for curved-wing — body junctures and finite wing tips. In 
addition, the right- and left-hand wing panels are considered separately. 

Pointed-tip wings .- Consider the right-wing panel shown in the fol- 
lowing sketch: 



(c) 
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The frontal projection of the area of one wing panel cut by the Mach 
plane is given by 



dy 


and - can be written as 
c 


z_ 

c 


z max z 
c z max 


1 1 


Then, for ^ constant, 

A - s o c oi / f(v > 

The value of tj and v are related by the intersection line of the Mach 
plane and the wing chord plane for the right-wing panel given by the 
equation 


x' = x Q + y(p cos 0 - tan A) 
and for the left-wing panel by the equation 

x' = Xq + y(p cos 0 + tan A) 


Then, 


With 


v = 


v~ + — (6 cos 0 ± tanA)t] 

0 c o 


£ - 1 - 1 


'] 


c 0 


v = 


— k + — 

- iL° c o 


(p cos 9 ± tan A)q 


1 


and 


Let 


= 


VQ - V 


(tan A ± p cos 0)— - v 
c 0 


s n 

m = tan A 



In terms of tapered wing geometry > m is given by 


m = A 1 + tan A 
k 1 - A 


Then, 


n = 


v 0 - v 


mil x *= — — -) - v 

tan A / 


,(i± 


P cos 6 ^ 


Let 


T , /■, . 3 cos 0 

K = m 1 + — 

\ tan A j 


for the left-wing panel and 


K = ml - 


3 cos 8 
tan A 


for the right-wing panel. Then, 



dr} K - v Q 

dv (K - v ) 2 


C _ 1 r, _ K - V 0 

55-1-1 -mr 

JL d, - av 

C ° (K - v) 5 

The equation for the area can then be written as 


A - -Vo | ( K - v o) 


2 P V upper f (v) 

J V. (K - v)3 

lower ' 


dv = 
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The slope of the area -distribution curve is obtained by differentiating 
the expression for A 


dA 

dv. 


'0 


" S 0 C 0 c 


-( K - »o) 


upper f( v ) 


lower 


(K - v) 3 


dv + 


(k - vq ) 2 

” V lower) 


dv 


lover 


dv, 


( v lower ) 


(K - v 0 ) 2 dv. 


upper 


or 


0 


dA 


'K - v y dv n 
l upper ) 0 


f ( v upper) 


dv 0 - s o c o I H(K jVt> ) 

Curves of G(K,vo) and H(K,vo) have been made up for a 65 A series 
airfoil and are given in figures l6 and 17 for values of K from 0 

to 2.4 and vq from 0 to 1. In evaluating f — — ■ - »» — dv, f(v) was 

^ (K - v) 3 

assumed to vary linearly between airfoil ordinate stations. Figure 18 
gives a plot of f(v) for this assumption. 

For K and Vq greater than 1, G(K, vq) and H(K, Vq) are given 
by the expressions 

G(K,v 0 ) = G(K,1) — " VQ) ' 

(K - l) 2 


H(K,vo) = H(K,1) 


(k - vq) 
(K - 1) 


For K and Vq less than 0 , G(K,vo) and H(K,vo) are given by 
the expressions 


G(K,vo) = G(K,0) ( - K " VQ ) 2 

¥L d 

h(k,v 0 ) = H(K,0) (K ~ VQ ) 

K 



23 


Correction for curved-wing — body juncture .- The following sketch 
shows a pointed-tip wing mounted on a curved body. 



(a) 

The areas and slopes will be referred to the actual wing geometry 

(c r , s, and >>). The areas and slopes, however, will be for the exposed 

pointed-wing tip. 

In sketch (d) consider one point of intersection of the wing panel 
with the body. The area of the wing cut by the Mach plane through this 
point is determined only by the product of s q c o the ex P ose ^ wing 

through the point and the value of v Q for the exposed wing at the point 
of intersection. As the point of intersection changes, s Q , cq, and v 0 

change and account for the intersection line. Expressed in terms of the 
actual wing-body characteristics, s q c q gtven by 


S 0 C 0 ” 


c r s 
1 - A 



2 


The quantity |j- 


is related to v Q by the expression 


v 


0 


v r - (1 - A)m J 
1 - (1 - A)§ 


The area of the exposed wing panel cut by the Mach plane can be written 
as 
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C r S “ r— -i 2 

A = — -\l - (1 - A){| G(K,vo) 

The area is calculated for given value of Vq. The center-line value of 
v is given by 

v r = v b[ 1 * d - *>§] + K(i - A)§ 

The slope is obtained by differentiating the expression for A and 
is given by 


c r s |[l - (I - A)f] 

[l - (1 - A)|]h(k,v 0 ) - 2(1 - A)^ G(K,vo) 

1 - A 4 

1 - (1 - A)§ + (K - v 0 )(l - A)!! 


dv 0 

“ rrlC 1 - (1 - 


Correction for finite wing tip .- In order to correct the pointed- 
tip wing panel and slopes for the finite wing tip, the areas and slopes 
outboard of the wing tip are subtracted. 



From sketch (e): 


Vo = 


c r s 


1 - A 



4X 


Then the areas and slopes are given by 


A. 


tip 


c r s £ A 2 
1 - A 


G(K,vo) 


^ip c r s | * 


dv 


h(k,vq) 


■ r 1 - A 
The center-line value of v is given by 
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v r = v q a + K(1 - A) 
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(b) 0 = 0°. 


Figure 1.- The areas and pressures which influence the drag of configura- 

tions at supersonic speeds. 
















Figure 6.- Comparison of the drag of parabolic bodies calculated with the area rule with 

experiment and slender-body theory. 












Tota/ 



;ure 7*“ Comparison of the area-distribution- curve slopes for a delta 
wing obtained by using the thin- and thick- wing assumptions. NACA 
65 AOO 6 airfoil section; M = 1. 4l4; 9 = 90°* 


L-1000 








/3cOi+/tanJL 


(a) Area distribution drag. 



(b) Area distribution drag against roll angle. 



(c) Configuration drag. 

Figure 10.- An example of the calculation of the configuration drag from 
the drag of the area distributions at various values of (3 cos 0/tan A 
60° delta wing: NACA 65 AOO 6 airfoil section. 
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AS A4 AS A6 


(a) Basic body; model 1 (ref. 16); 


fb 

Sw 


0.0305. 



'b) Model 4 (ref. 16); 6o° delta wing; NACA 65 AOO 3 airfoil section; 

Sb 


= 0.0305. 


w 



(c) Model 5 (ref. 16) ; 60° delta wing; NACA 65 AOO 6 airfoil section; 


<r- = 0.0305. 
^w 



(d) Model 6 (ref. 19) ; 60 ° delta wing. Thickness ratio varies from 0.03 


at root to 0.06 at 0.9 semispan. 


Sw 


0.0305. 


Figure 11 .- Comparison of the calculated drag with experiment and 
equivalent-body test results for wing-body combinations. 


onnT-l 







(h) Model 2 (ref. l6)j A = 3«04; A = 0.39^; ^c/k ~ ’ 'c ~ 

- 2 - = 0.0606. 

S w 

Figure 11.- Concluded. 



Figure 12.- Comparison of equivalent-body drag and configuration drag 

for an airplane configuration. 







45 



Figure 1 6 .- Area distribution parameter G(K, vq) for 65 A series airfoil. 


kb 



(a) vq from 0 to 0 .k^>. 

Figure 17.- Area-distribution- slope parameter H(K, vq) for 65A series 

airfoil. 






